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Abstract 

Observational results from several experiments like COBE, SP, Saskatoon, PYTHON, 
ARGO, MAX, MSAM, Tenerife and CAT are considered and a comparison is made with 
predictions from several models. Conclusions are reached about the viability of current 
structure formation models. 



1 Introduction 

Cosmic microwave background (CMB) radiation anisotropies provide a direct probe of the 
structure of the universe at early times and constrain the values of key cosmological parameters 
such as Qq, H q and Qf,. Experiments are now producing reliable measurements of the anisotropy 
on a variety of angular scales. They provide useful information about the process of structure 
formation, avoiding the systematic uncertainties and biases of conventional techniques. For 
Gaussian random fields all the statistical information is included in the angular power spectrum 
of the CMB temperature fluctuations. This being so, to determine the shape of this angular 
power spectrum constitutes one of the main goals of observational cosmology. We proceeded 
with a intercomparison of data and models in order to get the best-fit shapes of the angular 
power spectrum and constrain the model parameters. We present an extension of the work 
carried out by Hancock et al. [|12|, in which we used an analytic approximation to current 
CMB models and compared them with recent data. Here we consider some of the exact COBE 
normalised angular power spectrum models for the four-year COBE data, and apply the analysis 



to 4 of the 8 binned angular power spectrum points |26j for the COBE experiment in conjunction 



with Tenerife, Python, South-Pole, Saskatoon, MAX, ARGO and CAT experiments. We find 
that low values of Q are excluded by these exact models, in agreement with the results presented 
by Hancock et al. ||12|| . For a flat standard Cold Dark Matter dominated Universe we use our 



results in conjunction with Big Bang nucleosynthesis constraints to determine the value of the 
Hubble constant as Ho = 30 — 70 km s _1 Mpc~ x for baryon fractions Vtb = 0.05 to 0.2. For Ho = 
50kms _1 Mpc^ 1 we find the primordial spectral index of the fluctuations to be n = 1.1 ± 0.1, 
in agreement with the inflationary prediction of n ~ 1.0. Flat models with a contribution 
from a cosmological constant A seem to be consistent with the data. Models in which the 
fluctuations originate from cosmic strings appear to be consistent with the observations if we 
do not consider the higher calibration Saskatoon data. Textures model are only allowed for the 
lower calibration Saskatoon data. A more extensive discussion of these models will appear in 
Rocha et al. [[L9| (in preparation). 



2 Models 

Models of structure formation j|, [13| predict the shape and amplitude of the CMB power 
spectrum and its Fourier equivalent, the autocorrelation function C{6) = (AT(ni)AT(n 2 )) 
where n\.n2 = cos8. The intrinsic angular correlation function C(8) may be expanded in 
terms of spherical harmonics: C{9) = Ya>2^ + ^)CiPi(cos(6)) / Arr , where low order multipoles 
I correspond to large angular scales 9 and large Z-modes to small angular scales. The coeffi- 
cients Ci are predicted by the cosmological models and contain all of the relevant statistical 
information for models described by Gaussian random fields Q]. We consider the following 
models: flat CDM models with fi fc =0.01-0.3 and F =30-100 km s" 1 Mpc" 1 (G. Efstathiou, 
priv. comm.); open models with values for (f2o,^o) (where ho = ^/(lOOkms" 1 Mpc -1 )) of 
(0.1,0.75), (0.2,0.7), (0.3,0.65), (0.4,0.65), (0.5,0.6) (N. Sugiyama, priv. comm.); models with 
a cosmological constant, A and no reionization, with (A,h Q ,Qb) of (0.6,0.5,0.05), (0.6,0.5,0.03), 
(0.6,0.8,0.06), (0.7,0.5,0.05), (0.2,0.5,0.05); a string model (J. Magueijo, priv.comm.); and tex- 
ture models with (Q,h ,Q b ) = (l, 0.5, 0.05), (1,0.7,0.05) (N. Turok, priv. comm.). 

We also consider the case of tilted models with the spectral index n of the primordial 
fluctuations varying in the range n = 0.7-1.4 with (/i ,fi fe ) = (0.3,0.2), (0.45,0.1), (0.5,0.07), for 
direct comparison with scale invariant flat CDM models. 



3 Data 

We consider all the latest CMB measurements, including new results from COBE, Tenerife, 
MAX, Saskatoon and CAT, with the exception of the MSAM results and the MAX detection in 
the Mu Pegasi region ( for details on data selection see Hancock et al. [0). We also consider 
all 8 multipole bands of the angular power spectrum extracted from the 4-year COBE data 
PEfl . In order to avoid eventual correlations between these data points we have applied the 
fitting analysis to the 4 even or 4 odd multipole bands in conjunction with the remaining data 
points. The results are similar whether we use the even or the odd multipoles. The use of 4 
instead of one multipole is particularly useful in discriminating the value of the spectral index n. 
The different experiments sample different angular scales according to their window functions 
|28| , f27f , Wi, as shown in Fig. [I]. On the left hand side we present the window functions for 4 of 
the 8 multipoles of the COBE experiment. The conversion of data to a common framework in 



terms of flat bandpower estimates 0, pi follows the same procedure as in Hancock et al. [12]. In 



Table 1 we display these bandpower estimates AT; ± a where // and l u represent the lower and 



upper points at which the window function for each configuration reaches half of its peak value. 
For each experiment the uncertainties of the intrinsic anisotropy level have been calculated 
using a likelihood analysis (see e.g. Hancock et al. [ 10|| ) , which incorporates uncertainties due 



to random errors, theoretical uncertainty (cosmic variance) |2l|, ^] and enhancement of cosmic 
variance due to partial sky coverage (sample variance) ||23|| . The errors in AT; quoted in column 
3 are 68 % confidence limits and have been obtained by averaging the difference in the reported 
68 % upper and lower limits and the best fit AT;. This procedure introduces a small bias 
into the results since the likelihood function is in general only an approximation to a Gaussian 
distribution. 



4 Comparison 

4.1 Method 

For a given theoretical model it is possible to compute the value of AT b s one would expect 
to observe using the angular power spectrum CVs as predicted by the model and the window 
function corresponding to the experimental configuration Wi. We convert this to the bandpower 
equivalent result AT; and compare with AT° bs . We then compute the chi-squared for this set of 
C{s for the nd data points in Table 1, and the relative likelihood function is formed according 
to: L oc exp(— x 2 /2) ( for details see Hancock et al. [0]). As mentioned above the process by 
which the error bars quoted in Table 1 are computed introduces a small bias into the results. 
In order to assess the significance of this we applied the fitting analysis to an alternative set of 
data which consisted of the mean value of the data and corresponding error bars. We find no 
significant difference in the results obtained. We also vary the power spectrum normalization 
within the 95% limits for the COBE 4-year data [Q]. The data included in the fit are those from 
Table 1, which, with the exception of Saskatoon, include uncertainties in the overall calibration. 
There is a ±14% calibration error in the Saskatoon data, but since the Saskatoon points are not 
independent this will apply equally to all five points [p7| . The likelihood function is evaluated 



for three cases: (i) that the calibration is correct, (ii) the calibration is the lowest allowed value 
and (iii) the calibration is the maximum allowed value. We consider that a given model offers 
an acceptable chi-squared fit when P(x 2 ) > 0.05. 

4.2 Results 

We have applied the analysis to the COBE normalised flat CDM models (provided by G. 
Efstathiou) considering only one multipole for COBE namely C2 obtained from the value 
of the Qrms-ps- The models offering an acceptable chi-squared fit (P(x 2 ) > 0.05) to the 
CMB power spectrum whilst simultaneously satisfying nucleosynthesis constraints []5| , compre- 
hend 0.05 < £lb < 0.2, 30kms~ 1 Mpc _1 < H < 50kms _1 Mpc -1 . Considering the highest 
Saskatoon data calibration the constraints become 0.1 < < 0.2, 30kms _1 Mpc -1 < Ho < 
35kms _1 Mpc -1 . Allowing for the lowest Saskatoon data calibration relaxes the constraints up 
to H = TOkms^Mpc" 1 and 0.02 < Q b < 0.2 (see Hancock et al. [0). In general, recent 
optical and Sunyaev-Zel'dovich observations of the Hubble constant [IB, |8], imply Ho in 

the range 50 — 80kms _1 Mpc -1 . 

Applying this analysis to the open models we conclude that Q < 0.3 is not compatible with 
the data, while the value of Q = 0.4 is only allowed for the lower Saskatoon calibration data. A 
value of Q = 0.5 offers an acceptable chi-squared fitting for both case (i) and (ii) calibrations. 
This confirms the results obtained using the approximate formula of Hancock et al. reassuring 




us that the generalized parametrized formula previously used constitutes a good approximation 
given the uncertainties in the data. 

The models with low Q and a non-zero cosmological constant A seem to be consistent with 
the data, with large values of A providing the best fit. In Fig. |] left hand side we present the 
best candidates of both the open and A models compared with a flat CDM model with ho = 0.45 
and Qf, = 0.1 assuming case (i) calibration. Models in which structure formation is initiated 
by cosmic strings [16] are consistent with the data, but are excluded for the higher calibration 
Saskatoon data. Texture models are also only allowed if one considers the lower calibration for 
the Saskatoon data. In Fig . ||] right hand side we present a cosmic string model and the best 
candidate of the texture models (with Q = 1, h = 0.5 and fib = 0.05) for comparison with the 
flat CDM model of Fig. |2] using case (i) calibration. 

We also consider some of the best candidates of the flat CDM models: (1) h = 0.3, fib = 0.2, 
(2) h Q = 0.45, Qb = 0.1 and (3) h = 0.5, Qb = 0.07 in which we allow for the tilting of the power 
spectrum to vary within n = 0.7 — 1.4. The fitting for the spectral index n of the primordial 
fluctuations was done considering two sets of data: (a) one multipole, (b) 4 multipole bands 
for the COBE experiment and case (i) for the calibration. For case (a) we get for model (1) 
a value of n = 0.94 ± 0.06, for model (2) n = 1.03 ± 0.07 and for model (3) n = 1.08 ± 0.07. 
For case (b) we find for model (1) n = 0.99 ± 0.07, for model (2) n = 1.08 ± 0.07 and for 
model (3) n = 1.11 ± 0.07 (68 % confidence interval). Allowing for the three calibration cases 
and taking extreme limits the uncertainties are altered as follows: for case (a) for model (1) 
n = 0.94 ± 0.1, model (2) n = 1.03 ± 0.1 and model (3) n = 1.08±g;?| while for case (b) for 
model (1) n = 0.99tg;^, model (2) n = 1.08 ± 0.1 and model (3) n = l.ll±g;??. In particular 
fixing H = 50kms _1 Mpc _1 we find n = 1.1 ± 0.1 (68 % confidence interval) allowing for 
the three calibration cases and taking extreme limits. This tight limit rules out a significant 
contribution from gravity wave background for these models, in the case of power law inflation, 
but is consistent with the prediction of n ~ 1.0 for scalar fluctuations generated by inflation. 
In Fig. |3] left hand side we present the model (3) and its 68 % confidence interval assuming the 
case (i) for the calibration. On the right hand side is displayed the likelihood function for the 
parameter n for model (3) which peaks at n = 1.1 showing a significant estimate of the spectral 
index parameter based on the actual CMB data. 



Table 1: Detai 



s of data results usee. 



Experiment 



AT, (/iK) 



a (/iK) 



Ref 



C0BE2 

C0BE4 

C0BE6 

C0BE8 

Tenerife 

PYTHON 

South Pole 

ARGO 

MAX GUM 

MAX ID 

MAX SH 

MAX PH 

MAX HR 

Saskatoonl 

Saskatoon2 

Saskatoon3 

Saskatoon4 

Saskatoon5 

CAT1 

CAT2 



25.4 
28.1 
22.3 
31.9 
34.1 
57.2 
39.5 
39.1 
54.5 
46.3 
49.1 
51.8 
32.7 
49.0 
69.0 
85.0 
86.0 
69.0 
50.8 
49.0 



8.1 

3.9 

4.2 

24.3 

12.5 

16.4 

11.4 

8.7 

13.6 

17.7 

19.1 

15.0 

9.5 

6.5 

6.5 

8.9 

11.0 

23.5 

15.4 

16.9 



3 
5 

11 

25 

20 

91 

57 

95 

145 

145 

145 

145 

145 

86 

166 

236 

285 

348 

396 

608 



3 
5 

10 

22 

13 

50 

31 

52 

78 

78 

78 

78 

78 

53 

119 

190 

243 

304 

339 

546 



12 

8 

15 

30 

31 

107 

106 

176 

263 

263 

263 

263 

263 

132 

206 

274 

320 

401 

483 

722 
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Figure 2: The data points from table 1 compared to the exact forms of the C\ for on the left 
hand side: an Q — 1, ho — 0.45 and Qf, = 0.1 standard CDM model (bold line), an fl = 0.5, 
ho = 0.6 open model (dot line) and a flat £7a = 0.7, ho = 0.5 and Qf, = 0.05 model (dashed 
line). For on the right hand side: the same standard CDM model compared with a cosmic 
string model (dashed line) and an Q — 1, ho — 0.5 and fib = 0.05 texture model (dot line). 



<1 





Figure 3: Tilted model with h = 0.5, flj, = 0.07, n — 1.1 and its 68 % confidence interval 
(left hand side), likelihood function for the parameter n (right hand side). 
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